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Abstract 

We have obtained deep J^s-band images centered on QSO UM402 {zem = 2.856) 
using IRCS camera and A036 system at the Cassegrain focus of Subaru Telescope. 
A faint nearby galaxy (m/^ = 21.91 ± 0.26 in the Vega system) that lies ~ 2'.'4 north 
of the QSO sightline has been clearly resolved by this high resolution imaging. The 
non-detection of this object by the previous deep R-band observation indicates that 
it has a red color {R — K)vega > 3.3. This galaxy has irregular morphology with two 
close components (separation ~0''3). Given the small impact parameter {b= 19.6 kpc, 
at ziis = 2.531), it might be a candidate galaxy giving rise to the Lyman Limit system 
absorption at Zabs = 2.531 seen in the QSO spectrum. Careful subtraction of the 
PSF from the QSO image revealed the QSO host galaxy. We modelled in detail the 
host galaxy properties using the 2-D decomposition algorithm GALFIT, and found 
that QSO UM402 is hosted by a giant elliptical of nik = 20.17 ±0.07, Sersic index 
n = 5.94±0.70 and of a scalelength ~4Kpc, as bright as other resolved hosts of radio- 
loud QSOs (RLQ) at similar redshift, although UM402 itself is radio quiet (RQQ). 
Meanwhile, this QSO falls within the scale relation on the black hole mass vs. host 
luminosity of the local and other high-z samples, showing no significant excess of 
Mbh /logLy ratio at 2 ~ 3. 
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1. Introduction 

In the past decades, the study of quasar absorption line systems has proved to be a 
powerful method for galaxy formation studies. Lyman Limit systems as the sibling of quasar 
absorption line family, which are defined to be optically thick at the Lyman limit (A < 912A) and 
of a neutral hydrogen column density N{HI) > lO^^cm"^, attract much attention within these 
two years. This is because recent cosmological numerical simulations presented a " cold stream 
accretion" in early massive hot haloes as the main mode of galaxy formation, and Lyman Limit 
systems (especially those at the peak epoch of galaxy formation 2 ~ 3) are considered to be 
important candidates for the study of such accretion mode (Prochaska et al. 2009, Dekel et al. 
2009, Fumagalh et al. 2011). 

On the other hand, the studies on the high redshift quasar hosts have also received 
increasing attention lately, since they open an important avenue to study the assembly and 
evolution of massive galaxies, in particular in relation to the growth of their central black 
holes. Although the bulge components of local galaxies and QSOs were initially thought to 
have a near-linear scaling relation between their central black hole masses and the masses of 
their host spheroids, which motivated the theoretical perspective on the growth of the central 
black holes and the host galaxy formation, we now know that things are probably not that 
simple (Kormendy & Richstone 1995, Magorrian et al. 1998, Silk & Rees 1998, Wang & 
Biermann 1998, Marconi & Hunt 2003, Haring & Rix 2004). Graham (2012) and Graham & 
Scott (2013) have recently revealed that there is a quadratic Mbh — Msphcroid relation at low- 
masses (Mbh < 10® M0), which would have impacts on our better understanding of the galaxy 
formation scenario. However, an analysis of high redshift QSOs would be one key to a robust 
check of the theoretical models (Monaco et al. 2000, Wang et al. 2000, 2003, Bian et al. 2003, 
Peng et al. 2006, Schramm et al. 2008). 

Despite the intense activity in both the domains mentioned above, the seeing limitation 
of the ground-based observations affects significantly the detection ability of faint galaxies close 
to QSO sightlines, while the detection and analysis of QSO host galaxies at high-z is difficult 
even from space due to the sensitivity and the difficulty in PSF subtraction (Bahcall et al. 
1994, 1995, Mcleod & Rieke. 1995, Kukula et al. 2001). 

Thanks to the adaptive optics technique which offers a high spatial resolution for a 
powerful detection of faint companions or intervening galaxies of QSO nearby fields by con- 

* Based on data collected at the Subaru Telescope, which is operated by the National Astronomical 
Observatory of Japan. 
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centrating the flux in the core of the images. The abihty to determine the host morphology 
is much improved recently with the 8-lOm ground-based telescope equipped with AO systems, 
especially, if observing in the infrared which might minimize the difference in luminosity be- 
tween the host and nucleus again. Although correct decomposition of the compact nucleus and 
the extended host galaxy for high-z QSO remains a challenge, due to the difficulties in precise 
PSF determination, several ground- and space-based high resolution imaging of high-z QSOs 
have successfully resolved the host galaxies up to 2; > 2 (Aretxaga et al. 1998, Lacy et al. 2002, 
Groom et al. 2004, Hutchings et al. 2003, Jahnke et al. 2004, Kuhlbrodt et al. 2005, Falomo 
et al. 2008, Schramm et al. 2008). 

UM402 is a bright, high redshift radio-quiet QSO discovered by Macalpine & Lewis 
(1978). It is bright enough to permit detailed spectroscopic observations with a resolution 
ranging from a few to hundreds kms~^, and showing clearly the strong and sharp Lyman a and 
CIV emission lines, as well as a Lyman Limit system at 2 = 2.531 (Nhi > 4.6 x lO^^cm^) even 
with a low- resolution spectrum (Sargent et al. 1989). Previous deep imaging in the optical 
band of this QSO field reported the detection of several close neighbors {6 ~ 4''7 — 7"), which 
are all spectroscopically confirmed at redshift z <1 (Le Brun et al. 1993, Guillemin & Bergeron 
1997). Thus, the galaxy counterpart of the LLS seen in the QSO spectrum is still not identified, 
indicating it might be much fainter than the previous detection limit of mr{3asky) = 25.2, or 
much closer to the QSO sightline. 

UM402 again is one of the high-z QSOs which we selected for a pilot study on their 
host galaxies. Meanwhile, its image quality is the best for the observing run on Sept. 17-19, 
2003 (UT), and reasonably good for such a study. Several issues were carefully considered during 
the target selection, 1) high-z QSOs near the era of peak QSO activity and cosmic star formation 
history at 2; ~ 2 — 3, are specially selected due to their importance on the understanding of the 
galaxy formation scenario; 2) there should be a bright guide star {R< 15) sufficiently close to 
the QSO sightline (< 30"), in order to be observed with IRCS+A036 on Subaru telescope; 3) 
the emission lines Ha (6563 A ), H/3(486lA ), OII(3727A ) and OIII(5007A ) should be avoided 
to be included in the observing bands. This is important for the host mass estimation and the 
host continuum property studies; 4) there exists a suitable PSF calibration star for the QSO 
host fitting in the further analysis. We will elaborate on this point in section 2.1. 

In this paper, we present the initial results from an adaptive optics(AO) deep imaging 
in the Ks bands of QSO UM402 at 2 = 2.856 with the IRCS camera on Subaru telescope. The 
cosmological parameters Q = 0.27, A = 0.73 and H = 71km/s/Mpc are adopted throughout. 
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2. Observations and data reduction 

2.1. Selection of AO guide star and PSF calibration star 

For Subaru A036 adaptive optics systems, a natural guide star(NGS) sufficiently close 
to the siglitline of the target is required as a reference source to assess the degradation of the 
wavefronts due to the turbulent atmosphere. There is a bright star with magnitude of R = 13.8 
and an angle distance to UM402 of 9gs ~ 31". We selected this bright star as our AO guide 
star, and expected to obtain the AO corrected PSF better than 0'.'2 if the natural seeing is 
< 0% (Takami et al. 2004). 

Usually, the guide star could not be used to directly model the PSF. This is because 
the PSF is expected to change with the angular distance from the AO guide star, and the 
actual PSF at the position of the target will be degraded. On the other hand, there will be 
the saturation problem for a very bright guide star required by the optimal AO correction. In 
this case, we have to select other PSF calibration star which could be observed at a condition 
as similar as possible to that of the QSO, i.e. similar magnitude, similar direction and angular 
distance to the guide star. We selected a suitable PSF calibration star for UM402, which is 
of similar brightness, similar guide star distance, but a guide star angle offset ~ 180° to that 
of QSO (see Fig. 1 left). The details of the QSO, the guide star and the PSF reference star 
are listed in Tab. 1. The guide star used for AO correction is the same for QSO and the PSF 
calibration star. 

2.2. Observations and data reduction 

The AO-assisted Ks band deep imaging of UM402 was made on Sept. 17-19, 2003 (UT), 
using the IRCS camera on Subaru 8.2m telescope at Mauna Kea and the Subaru Cassegrain AO 
system with a 36 element curvature wave front sensor, as well as a bimorph-type deformable 
mirror with the same number of elements (A036; Takami et al. 2004). The camera uses one 
1024 X 1024 InSb Alladin III detector and has two imaging modes with different pixel scales. 
We adopted in the observation a pixel scale of 0''023 (23mas mode), providing a field size of 
23" X 23"(IRCS; Kobayashi et al. 2000). In order to remove the bad pixels, we adopted nine- 
point dithering in a 3 x 3 grid with a dithering step of 5", 6"5 or 7". To reduce the readout 
noise, 16 times nondestructive readout (16-NDR) was applied for each readout. Dark frames 
and dome fiats were taken at the end of each nights. Most of observing nights were clear and 
photometric. The median seeing size was ~ 0"5, and the airmass was mostly smaller than 1.4. 
After an optimal function of the AO system was achieved, we offset the FOV of the telescope 
to put the QSO or the PSF calibrator star in the center of the FOV. 

Similar as other currently available AO systems, we suffered from the small field view 
of the AO detectors, and were not able to include simultaneously a suitable PSF calibration 
star in the QSO exposures to directly evaluate the PSF. In order to monitor and assess the 
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Table 1. Observed QSO, PSF star and the guide star 

Obj. Type RA(J2000) DEC(J2000) z R^ag iexp(lirs) FWRM"^ GS(d)^ 

UM 402 RQQ 02 09 50.71 -00 05 06.6 2.855 15.8 4.52 0.13 14.54 31 

PSF star 02 09 54.51 -00 05 34.0 16.6 2.34 0.11 15.47 30 

Guide star 02 09 52.84 -00 05 15.2 13.8 

° Image quality measured as the FWHM for the coadded Ks images of ah good exposures in arcseconds. 

Observed K magnitude of the target for the coadded Ks images of ah exposures. 

^ Distance in aresec from the QSO (PSF star) to the Guide star. 

temporal variability of the PSF, we observed the PSF calibration star just before and after the 
observations of the QSO. More specifically, we observed the QSO itself in an exposure set of 
9 X 80s, or 9 X 70s using the 3x3 dithering pattern, nested between similar dithering observa- 
tions of the PSF calibration star. Such interleaving observations could provide us information 
on the temporal variation of the PSF during the target observations. Meanwhile, such non- 
simultaneous PSF calibration would provide very similar correction quality of the AO systems 
to the QSO images, since the variability applies to both the PSF star exposures and the QSO 
images displaying similar Strehl values. 

We adopt the "core width r2o", which is defined to encircle 20% of the total flux of 
a point source, as the image quality indicator based on a close relation between and the 
Strehl ratio given by Kuhlbrodt et al. (2005). We measured the value of the QSO and the 
PSF calibration star for each exposure frame of all three observing nights, and plotted out the 
progression of the "core width r2o" in case of 1) good exposures with r2o < 3.5 pix (Strehl ratio 
S > 50%, see Fig. 2 top); 2) good and stable images with r2o < 3 pix, as well as the first two 
dithering blocks on Sept. 18 were excluded (Fig. 2 bottom). Our further discussion on the 
host galaxy properties mostly relies on these good exposures. 

The package "IRCS-IMGRED" was used to make the dark frame, flat frame, bad pixel 
mask and sky frame, as well as fiat fielding and sky subtraction (Minowa et al. 2006). Finally, 
the dithered frames of the QSO and the PSF calibration star were aligned and averaged respec- 
tively, using an outlier rejection algorithm. The ^Tg-band image quality measured as the full 
width at half maximum (FWHM) for both the QSO and the PSF calibration star is presented 
in Tab. 1. 

The standard star p533-d and p338-c were observed as the photometric calibrator, which 
were selected from Hawarden et al. (2001). All magnitudes subsequently quoted in this paper 
are in the Vega system. 

2.3. PSF construction and simple subtraction 

Our goal with the high resolution AO images is to detect any possible faint galaxies 
nearby the QSO sightline, as well as the faint extended host galaxy hidden in the glare of 
the central bright QSO light. Therefore, it is mandatory to estimate properly the AO PSF 
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and subtract the light contribution from the bright central point source of the QSO images, 
in order to unveil the underneath faint objects, and to reduce the effects to their photometric 
measurements. 

We created a PSF using interleaved exposures of the PSF calibration star between QSO 
observations. A synthetic PSF was constructed by fitting the coadded PSF calibrator star 
image with a number of Sersic components of slightly shiftted peak position, magnitude and 
position angle. We kept increasing the number of Sersic components of the PSF star until the 
residuals for the star and the neighboring contamination (if any) go down to minimum, which 
means there is no obvious structural component left in the residual image, as well as not much 
improvement evenif we add extra Sersic models to the fitting. The Sersic components which 
were used to fit the PSF calibrator star were summed up into a single PSF image which was 
then carefully trimmed to have the PSF centered in the PSF image. More specifically, if the 
number of rows/columns along one side of the PSF star image (NPIX) is an odd number, then 
the PSF peak will be centered at NPIX/2 + 0.5. Otherwise, if NPIX=even, the PSF peak 
should be centered at NPIX/2 + 1. This synthetic PSF will be used later on for the 2-D image 
decomposition and would be very helpful for the GALFIT code to converge properly. 

To detect the host galaxy, we first subtracted the PSF from the QSO image using a 
very conservative and simple method, same as other studies on the high-z QSO host galaxies. 
We scaled the PSF flux to the QSO central peak intensity and aligned them. In this way, 
such a method implies an oversubtraction of the nuclear component from the inner region, and 
provides a model independent host detection and a low limit on the host flux (Sanchez et al. 
2004). The contour plots of the PSF subtracted good QSO images for each observing night 
from Sept. 17 - Sept. 19 were given in Fig. 3 (left) , while Fig. 4 (left) shows the contour plot 
of the PSF subtracted coadded good QSO images of all three nights. Further one-dimensional 
radial profile analysis using ELLIPSE within IRAF, as well as two-dimensional model fitting 
using GALFIT on the host properties will be presented in section 3.2. 

3. Analysis 

3.1. Galaxies near the QSO sightline and the luminosity 

The high resolution deep image has clearly resolved two galaxies within a 5"x 5" field 
around UM402 (see Fig. 1 right). They are a fuzzy galaxy about 4''7 southern of the QSO, 
and a close object ~ 2''4 northern of the QSO sightline which appears as a double system 
with a separation of the two components ~ 0''3. For the left component of this double system, 
there seems to be a faint tidal-tail like feature towards southeast, suggesting a possible merging 
system for this object. The fuzz galaxy southern of the QSO sightline has been detected by 
SDSS and other optical deep imaging as a nearby irregular galaxy at 2; ~ 0.36 (Le Brun et al. 
1993, Guillenmin & Bergeron 1997). 
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We measured the photometry of this close double system by running a 2-D decomposition 
algorithm GALFIT on the QSO image, where the QSO(psf), the host galaxy and any other 
nearby galaxies (Sersic profiles) in the field were fitted simultaneously to deblend everything 
together, in order to reduce the contaminating flux from the wing in the PSF of the QSO 
(Peng et al. 2002). The best fitting gives an apparent magnitude = 21.91 ±0.26 for the 
neighboring galaxy which lies ~ 2'.'4 northern of the QSO sightline. Guillenmin & Bergeron 
(1997) did not see this northern object, though their image quality would allow to resolve it 
from the QSO light. This might mean that it is fainter than the detection limit of the previous 
imaging observations (m^ = 25.2), and of a red color [R — K)yega > 3.3. 

According to the impact parameter vs column density relation {b — logNni) for all 
confirmed DLA and LLS absorbers given by MoUer & Wallen (1998), we suspect this close 
double system might be a candidate galaxy giving rise to the Lyman Limit absorption at 
Zabs ~ 2.5 previously seen in the QSO spectrum {Nhi > 4.6 x 10^ ''cm^). Comparing to the 
apparent K-band magnitude vs. stellar mass relation for objects at 2.3 < z < 2.6 from the 
MOIRCS Deep Survey (MODS), this galaxy would have a stellar mass comparable to the 
median stellar mass (~ a few 10^° Mq) of the MODS sample (Tanaka et al. 2011). 

Further observations for the spectroscopic redshift of both components of the double 
system are strongly required. If confirmed, it would be an important high-z evidence of a 
merging system as Lyman Limit absorber. 

3.2. Host galaxy detection and the luminosity 

We measured the radial profiles of the coadded good images of the QSO and the PSF 
star (r2o < 3.5 pix) using STSDAS task ELLIPSE, after masking out the close companions. The 
radial profiles of the coadded good QSO images of each observing night are compared with 
the PSF stellar profiles acquired during QSO exposures of that night (Fig. 3 right), while the 
radial profiles of the coadded good QSO and PSF exposures of all three nights, as well as their 
residual are presented in Fig. 4 (right). The extended emission of the QSO host galaxy at the 
radii > 0''2 is resolved for almost all three nights, although the images obtained on Sept. 18 
are not good enough. 

We further modelled the 2D luminosity distribution of the host galaxy by fitting simul- 
taneously the QSO nucleus and its host galaxy, as well as any neighboring companions using 
GALFIT. With the constructed synthetic PSF image, we made a model fitting to the coadded 
good QSO image, where the QSO nucleus (psf), the host galaxy and the neighboring galaxies 
(Sersic profiles) are fitted simultaneously to deblend everything together. The best fit results 
for the coadded good images of all three observing nights (x^ ~ 1) indicate that the host galaxy 
is a luminous elliptical of = 20. 17 ±0.07, Sersic index n = 5.94 ±0.70 and the effective radius 
Re = 0'.'58 ± O'.'OQ (case c in Tab. 2, or the last line of Tab. 3). 

We compared the host luminosity of UM402 with other high-z QSO hosts at similar red- 
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shift, and found that it is well located in the range of the observed apparent K band magnitude 
versus redshit (the K-z plot) for luminous radio galaxies and other radio-loud QSO hosts of 
similar redshift range, although it is radio quiet (see Fig. 5). Our detection of the QSO host 
galaxy in this analysis demonstrates that radio-quiet QSO host of bright elliptical morphology 
exists up to 2; ~ 3, similar as what we see at low redshift (McLure et al. 1999). 

Since at 2; ~ 3, the observed i^^-band corresponds approximately to the rest-frame V- 
band, we can make an easy comparison of our QSO host galaxy {My = —25.28 ±0.07) with the 
average radio-quiet QSOs at z < 0.46 (< My >~ —22.7) (Smith et al. 1986, Hutchings et al. 
1989, Bahcall et al. 1997, Hamilton et al. 2002). Thus, the host of UM402 is likely to be about 
3.0 mag brighter than the low-z samples. Again, if we consider that a local elliptical has 
a luminosity of My = -22.35 (Hamilton et al. 2002), The host galaxy of UM402 is about 3.0 
mag brighter than a non-evolving elliptical, but not significantly brighter than a passively 
evolved one if placed at 2 ~ 3 (Barger et al. 1998, Bruzual & Chariot 2003, Miyazaki et al. 
2003 ). 

From our best 2D fitting, the effective radius of the host galaxy of UM402 is ~ 4 Kpc 
[Re = 0''58±0''06 at Zem = 2.856), as well as a stellar mass ~ 10^^ M0 based on the relation 
between the stellar mass and the observed total Ks magnitude for galaxies at 2 < 2; < 3 in the 
FIRES, GOODS and MUSYC fields (van Dokkum et al. 2006). Considering that SDSS galaxies 
on the red sequence and of similar stellar masses (~ 10^^ Mq) would have a median scalelength 
-Re ~ 5 Kpc, the host galaxy of UM402 is only slightly smaller than the local massive galaxies, 
indicating that a massive host galaxy of elliptical morphology could be already formed as early 
as 2 ~ 3. However, it is clearly larger than the median scalelength of a sample of quiescent 
massive galaxies at similar redshift presented by Damjanov et al. (2011). 

3.3. PSF variability and the systematic errors 

PSF uncertainty is always a stressful concern of the model fitting on the galaxy profiles, 
especially for the compact sources, which requires very good knowledge of the PSF. To check the 
systematic errors due to the PSF variation, we performed four different empirical approaches 
as following: 

1) we masked out the central regions of the QSO images which cannot be fitted well, and did 
again similar PSF fitting to the coadded QSO image. No systematic errors were found. 

2) we grouped the quasar and the PSF star exposures into three different sections according 
to their observing condition, i.e. a) very good and stable images with < 3 pix, as well 
as the first two dithering blocks on Sept. 18 were excluded (see Fig. 2, bottom); b) very 
good exposures with < 3 pix; c) good exposures with < 3.5 pix, shown in Fig. 2 
(top). These three cases are dubbed as case a), case b) and case c) in Tab. 2. We 
performed similar PSF fitting to the coadded images of above three cases respectively, and 
presented the fitting results in Tab. 2. 
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Table 2. 2-D decomposition results for coadded QSO image of different observing conditions 

Obs. Date PSF niag.(mfc) host mag.(mfe) effective radius (") Sersic index b/a 

case a) 19.56±0.02 20.20±0.22 0.66±0.40 *20.0±6.9* 0.94±0.04 1.07 

caseb) 19.62 ±0.01 20.17 ±0.09 0.58 ±0.08 5.94 ±0.86 0.85 ±0.02 1.3 

casec) 19.62±0.01 20.17±0.07 0.58±0.06 5.94±0.70 0.85±0.02 1.24 

n3-all 19.62 ±0.05 20.17 ±0.08 0.58 ±0.06 5.94 ±0.88 0.85 ±0.01 1.08 



3) a similar GALFIT model fitting as in approach 2) was applied to each observing night 
respectively. The fitting results are shown in Tab. 3. We listed at the last line of Tab. 3 
the results of fitting on the coadded good images of all three nights (i.e. case c), to show 
how S/N ratio would affect the fitting results. 

4) to see how sensitive the GALFIT results can be to the adopted PSF, we ran GALFIT 
model fitting using a range of PSFs to assess the uncertainty associated with the PSF. 
The set of PSFs were constructed using each dithering blocks (9-dithered frames) of the 
PSF star observed in good and stable conditions (i.e. case a in Tab. 2, or Fig. 2, bottom). 
We named the set of PSFs according to their dithering sequence and the observing date, 
such as, nl7-ld and nl7-2d are the first and second dithering blocks of the PSF star observed 
on Sept. 17. We showed the radial profiles of the coadded good and stable QSO images of 
each observing night (red plus), compared with the PSF stellar profiles (green cross) acquired 
during QSO exposures and in imaging stacks of 9-dithered exposures in Fig. 6, as well as 
presented in Tab. 4 the results of model fitting on the coadded good QSO images (r2o < 3pix) 
of all three nights, using the set of good PSFs mentioned above. We found no significant 
systematics from such a trial. 

We noticed that the Sersic index from the GALFIT fitting on the combined good images 
of a single observing night on Sept. 17 and 18, as well as case a), goes unacceptably large 
(marked by * in tables). This is because GALFIT fitting is very sensitive to the sky background. 
A small error in the background would drive the index up to an unreasonably large number, 
leading to large systematic errors. However, the GALFIT fitting parameters for the observing 
night of Sept. 19, as well as case b) and case c) are reasonable and consistent, due to the better 
signal-to-noise ratio of the coadded good images. For a further check, we listed at the last line 
of Tab. 2, the GALFIT fitting on the combined images of all three nights, including those 
exposures with 3.5pix < < 4.0pix to improve the S/N ratio. We do see a reduced x^, while 
all other fitting parameters are still consistent with that of caseb) and casec). Considering 
all these empirical approaches, we think that the resolved host galaxy and the fitting results are 
reasonable, and were negaligibly affected by the temporal PSF variation during the observing 
run. 
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Table 3. 2-D decomposition results for coaddcd QSO image of each observing night and nights 



Obs. Date 


PSF niag.(TOfc) 


host niag.(mfc) 


effective radius (") 


Sersic index 


b/a 




Sept.17 


19.75 ±0.01 


20.38 ±0.25 


1.68±1.39 


*20.0±5.2* 


0.85 ±0.03 


0.95 


Sept. 18 


19.69 ±0.07 


20.60 ±1.95 


1.68±9.17 


*20.0±53* 


0.99±0.17 


2.2 


Sept. 19 


19.62 ±0.02 


20.17±0.26 


0.58 ±0.25 


5.94±2.51 


0.85 ±0.08 


8.11 


Sept. 17-19 


19.62 ±0.01 


20.17 ±0.07 


0.58 ±0.06 


5.94 ±0.70 


0.85 ±0.02 


1.24 




Table 4. GALFIT results 


using the set of good PSF 


s during this run 






Obs. Date 


PSF niag.(mfc) 


host mag.(mfe) 


effective radius (") 


Sersic index 


b/a 




nl7-ld 


19.62 ±0.01 


20.17 ±0.08 


0.58 ±0.07 


5.94 ±0.74 


0.85 ±0.02 


1.30 


nl7-2d 


19.57±0.01 


20.42 ±0.15 


1.20 ±0.35 


7.59 ±1.39 


0.90 ±0.03 


1.11 


nl8-3d 


19.62 ±0.01 


20.17 ±0.09 


0.58 ±0.08 


5.94 ±0.90 


0.85 ±0.02 


1.97 


nl9-ld 


19.62 ±0.01 


20.17 ±0.08 


0.58 ±0.08 


5.94 ±0.79 


0.85 ±0.02 


1.75 


nl9-2d 


19.62 ±0.01 


20.17 ±0.08 


0.58 ±0.08 


5.94 ±0.84 


0.85 ±0.02 


1.31 



3.4- Virial BH mass and the Mbh — Ly relation 

We estimated the black hole mass of QSO UM402 using the CIV broad emission line from 
its SDSS spectra. Based on the assumption that broad line region (BLR) is in approximately 
virial equilibrium and the scaling relation between luminosity and BLR size (Vestergaard 2002, 
McLure fcJarvis 2002, Kaspi et al. 2000, Mathur & Grupe 2005), Vestergaard & Peterson 
(2006) derived a prescription for the black hole mass measurement using broad emission line 
CIV as following: 

0.53' 



Mbh = 5.4 X 10^ 



FWHM 

km/s 



ALa(1350A) 



(1) 



lO^^erg/s 

where L;j(1350A) is the monochromatic luminosity at 1350 A, and FWHM is the full width at 
half maximam of the CIVA1550 emission line. 

We fitted the continuum using a power law in the emission-line-free windows: 1440- 
1460A, 1680-1700A, 2160-2180A. Since the red wing of CIV is contaminated by He II, only 
the spectral region of 1450 — I6OOA is considered in fitting the CIV line. After subtracting 
the fitted power law, three gaussian components were adopted to fit the selected range using 
the IDL mpfit package. The absorption lines were masked out during the fit. The best fit 
resuhs (x^ ~ 1.2) are FWHM=5457.13±345.39 km s"^ for the CIV line, and the flux density 
f^350=62.13(±0.32)xl0"^^ ergs s-^cm^^^-i. The black hole mass of UM402 is estimated to 
be logMbh ~ 9.9(±0.4dex) M0 (the SDSS systematic calibration uncertainty 0.36 dex has been 
included). It is a very massive black hole as we expected given its high luminosity. 

Shen et al. (2008) compiled black hole masses for ~ 60,000 SDSS quasars in the red- 
shift range of 0.1 ^ 2 ^ 4.5, where UM402 was included in their QSO list. However, only 
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FWHM=5628 km s for CIV line was given by them, which is consistent with our analysis 
here within the errors. 

In Fig. 7, we compared the black hole mass and the host luminosity of UM402 with those 
of other high redshift QSOs where their host galaxies were resolved and the black hole masses 
were estimated, to investigate where our object would fall within the Mbh — Ly relation. We 
found that the ratio of the black hole mass to its host luminosity of UM402 is in good agreement 
with the Mbh — logLy correlation of the low-z bulges and ellipticals given by Dunlop et al. 
(2003), and other high-z QSO samples (Kukula et al. 2001, Ridgway et al. 2002, Kuhlbrodt et 
al. 2005, Peng et al. 2006, Schramm et al. 2008). 

4. Summary 

We have presented an analysis of adaptive optics deep images in the Ks bands centered 
on QSO UM402 at Zem = 2.856, with IRCS camera and A036 systems on Subaru telescope. 

A close faint galaxy (~ 2''4 north of the QSO sightline) has been clearly resolved by 
the high resolution imaging, and appears as a double system with a separation of the two 
components ~ 0''3. The faint tidal-tail like feature from the left component of the double 
sytem indicates that it is probably a merging system. According to the empirical relation on 
the impact parameter vs. neutral hydrogen column density of all confirmed DLA and LLS 
absorbers given by Moller & Wallen (1998), as well as its red color (R — K)vcga > 3.3, we suspect 
that this faint object might be a candidate galaxy giving rise to the Lyman Limit absorption 
at Zabs = 2.531 previously seen in the QSO spectrum. If the redshift of both components of the 
double system are spectroscopically confirmed in the future, it would be an important high-z 
evidence of a merging system as the Lyman Limit absorber. 

After carefully subtracting the PSF, we are able to see the extended emission from the 
quasar host galaxy. We have run a 2-D imaging decomposition on the QSO nearby field of 
the coadded images with good quality (r2o < 3.5 pix) by GALFIT. The best-fit results indicate 
that QSO UM402 is hosted by a giant elliptical of rrik = 20. 17 ±0.07 and the effective radius 
of Re = 0''58±0''06. The host luminosity of UM402 is consistent with other resolved hosts of 
radio- loud QSOs at similar redshift range (Falomo et al. 2008). Since UM402 is a radio quiet 
QSO, this again indicates that radio quiet QSO host of bright elliptical morphology exists up 
to z ~ 3, same as what we see at low-z (McLure et al. 1999). However, it is in disagreement 
with the argument by Falomo et al. (2004) from a sample of 1 < 2 < 2, that RLQ hosts are on 
average a factor of ~ 2 luminous (massive) than RQQ hosts, a difference that does not appear 
to change significantly with the redshift. 

Further comparison of our QSO host with the average characteristics of radio-quiet QSOs 
at low-z and the local L^, ellipticals, shows that the host galaxy of UM402 is ~ 3 mag brighter 
than the low-z sample, not surprisingly brighter if considering a passive evolution of the host 
galaxy. Meanwhile, the scalelength of this QSO host {^AKpc) is only slightly smaller than 
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Fig. 1. Left: Finding cliart of QSO, guide star (GS) and PSF ealibration star: Riglit: A combined image 
of a region of ^ 5" x 5" around QSO UM402 (central briglit source) in the Ks band. The LLS absorber 
candidate is ~ 2'.'4 north of the QSO sightline. The pixel scale is 0'.'023 and the AO-corrected FWHM 
O'.'IS. North is up and East to the left. 

that of the local massive galaxies of similar stellar mass ~ 10^^ M0, but significantly larger than 
the median size of massive quiescent galaxies at similar redshift and of similar stellar mass 
(Damjanov et al. 2011). Although we need more samples to get better statistics, our results 
here might imply that QSO host galaxies with elliptical morphology would follow a different 
size evolution trend from normal early type galaxies. 

To see where our QSO falls within the Mbh — Ly relation, we estimated the black 
hole mass of this QSO using the CIV broad emission line from its SDSS spectra, adopting the 
prescription given by Vestergaard & Peterson (2006). We find that UM402 is in good agreement 
with the Mbh — logLy relation of the low-z bulges and ellipticals given by Dunlop et al. (2003), 
and other high-z QSO samples. 

Further information on the star formation history of a large sample of high-z QSO 
hosts, would help us to understand how AGN activity would affect the host galaxy formation 
and evolution, and make a comprehensive picture for the growth of central black holes and the 
formation of the host galaxies. 
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Fig. 2. Top: The "core width r2o" measured for each exposure frame of the QSO (green cross) and the 
PSF cahbrator star (red plus) vs. observing time sequence. We plot out here only good exposures with 
r2o < 3.5 pix for both the QSO and the PSF star; Bottom: Same as the top figure, but for good and stable 
exposures of both the QSO and the PSF star with r2o < 3 pix, as well as the first two unstable dithering 
blocks on Sept. 18 were excluded (indicated by the red arrowheads in the top plot). 
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Fig. 3. From top to bottom, we show here the contour plot of the extended emission of the quasar host 
revealed after PSF subtraction (left), as well as the observed radial surface brightness profile (right) of 
the QSO (red plus) and the PSF star (green cross), for the coadded good images (r2o < 3.5 pix) of each 
observing night from Sept. 17 to Sept. 19, shown in Fig. 2 (top). The contour plots shown in the left 
have a scale in pixels (O'.'023/pix). 
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Fig. 4. The contour plot of the extended host emission revealed after PSF subtraction (left), as well as 
the observed radial surface brightness profile (right) of the QSO (red plus) and the PSF star (green cross), 
for the coadded good images of all three nights. The residual of the two profiles (A/j,) is plotted with red 
plus at the bottom panel. The extended emission of the quasar host is resolved at the radius > 0'.'2. 
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Fig. 5. Apparent K band magnitude versus redshift (the K-z plot) for luminous radio galaxies (asterisks) 
(Willott et al. 2003), hosts of RLQ at z > 1 (circles for VLT data by Falomo et al. 2004 and Kotilainen 
et al. 2007; squares for HST data by Kukula et al. 2001 and open stars by Peng et al. 2006; pentagons 
for two RLQ by Falomo et al. 2008). The host of RQQ UM402 of this work is presented in the figure by 
a red filled star. 
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Fig. 6. Top: The radial profiles of the coadded good and stable QSO images (r2o < 3 pix) observed on 
Sept. 17 (red plus), compared with the PSF stellar profiles (green cross) acquired during QSO exposures 
and in imaging stacks of 9-dithered exposures (nl7-ld and nl7-2d); Middle: Same as the top figures, but 
for the good QSO and PSF star exposures obtained on Sept. 18; Bottom: Same as the top figures, but 
observed on Sept. 19. 
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Fig. 7. Relation between black hole masses and host luminosities. The red diamond is the result of this 
work. Other data points are taken from literature. The dotted line shows the linear relation of local 
sample from Dunlop et al. (2003). 
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